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A dark energy component is responsible for the present stage of acceleration of our universe. If no 
fine tuning is assumed on the dark energy potential then it will end up dominating the universe at 
late times and the universe will not stop this stage of acceleration. On the other hand, the equation 
of state of dark energy seems to be smaller than -1 as suggested by the cosmological data. We take 
this as an indication that dark energy does indeed interact with another fluid (we consider fermion 
fields ) and we determine the interaction through the cosmological data and extrapolate it into the 
future. We study the conditions under which a dark energy can dilute faster or decay into the 
fermion fields. We show that it is possible to live now in an accelerating epoch dominated by the 
dark energy and without introducing any fine tuning parameters the dark energy can either dilute 
faster or decaying into fermions in the future. The acceleration of the universe will then cease. 



I. INTRODUCTION 

A dark energy component is probably responsible for 
the present stage of acceleration of our universe [1], [2]. 
Perhaps the most appealing candidate for dark energy 
is that of a scalar field, quintessence [3], which can be 
either a fundamental particle or a composite particle [4] . 
Since dark energy dilutes slower than matter we expect 
it to dominate the universe at late times if no fine tun- 
ing is assumed on the dark energy potential V and the 
universe will not stop the present stage of acceleration. 
However, this conclusion can be overcome by assuming a 
dark energy interaction [5] . 

On the other hand, the equation of state of dark en- 
ergy seems to be smaller than minus one as suggested 
by the cosmological data [1],[2]. In general fluids with 
w < — 1 give many theoretically problems such as stabil- 
ity issues or wrong kinetic terms as phantom fields [11]. 
However, interacting dark energy [6]- [9], [5] is a very sim- 
ple and attractive option which we will use in this letter. 
The interaction must be quite weak since dark energy 
particles have not been produced in the accelerator and 
because the dark energy has not decayed into lighter (e.g. 
massless) fields such as the photon. 

It was suggested to determine the interaction between 
dark energy and this other fluid through the cosmological 
observations and extrapolate it into the future. Doing so, 
it was shown within a general framework that the dark 
energy dilutes exponentially fast and matter prevails in 
the end [5]. 

In this letter we take the dark energy to interact with 
fermion fields. It is tempting to take the fermions as 
neutrinos [8] since the order of magnitude of the energy 
of dark energy and neutrinos is similar. Furthermore, if 
the mass of neutrinos is larger than 0.8eV, as implied by 
the Heidelberg-Moscow experiment, then the dark energy 
cannot be a cosmological constant [9]. However, in this 
letter we do not assume the fermions to be neutrinos. 

We take the interaction between dark energy and 
fermions through a field dependent mass M = /(</>), 



where M is the mass of the fermions [8]. The function 
/(</>) plays an important roll in the evolution of the dark 
energy. The total potential Vr of the scalar field has a 
contribution of the original scalar potential V(cj>) and the 
interaction term f ((/))■ We calculate the conditions for ac- 
celeration and we show that the naive slow roll conditions 
Vt/Vt <C 1 and ViJ /V <C 1 do not imply acceleration. 
We determine the requirements under which a dark en- 
ergy can dilute faster or decay into the fermion fields. 
We obtain that only in the case where the mass of the 
dark energy, given by m 2 = V" + p-^f'/f, is dominated 
by the interaction term can the dark energy decay. We 
also show that it is possible to have dark energy diluting 
faster than the fermion fields. 

This letter is organized as follows. In section II we 
present our assumptions about dark energy in the ab- 
sence of interactions while in section III we present the 
dynamical equations of dark energy coupled to fermions 
via a fields dependent mass. In section IV we present the 
generic requirements for a particle decay and in section 
V we set the conditions for dark energy to decay in the 
future. In section VI we briefly show the generic evolu- 
tion of two coupled fluids and we discuss the evidence 
for dark energy decay. Finally, in section VII we study 
the possibility of dark energy diluting or decaying into 
fermions and we then present our conclusions in VIII. 



II. DARK ENERGY 

We will consider in this letter that the dark energy is 
given in terms of a scalar field with lagrangian 



(1) 



where V(<f>) is only a function of </>. We define the energy 
density and pressure of the scalar field as 



P<t> = \<t> 2 + V (</>), 



P4> = \ft ~ V(4>), (2) 
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i.e. without any interaction with other fields, and 

w _ P0 = ft/2 - V(<f>) 

W<P P4- ft/2 + V{4>) ' 

Clearly w<j, > — 1 at all times if > 0. For V to inflate 
the slow roll constrains to the potential must be satisfied, 



(3) 



\V'/V\ « 1, 



\V"/V\ « 1, 



(4) 



where a prime denotes derivative w.r.t. <fi, i.e. V = 
dV/d(f>. We do not want to introduce any arbitrary dark 
energy scale so we will assume that the slow roll condi- 
tions are satisfied also in the future and therefore V is a 
runaway potential and tends to zero at <f> — > oo. From 
eqs.(4) we see that V and V" also approach zero at late 
times and V is therefore negative. Inflation occurs in 
general for <j> > 1 and the mass m ~ H. 



III. INTERACTING DARK ENERGY 

The interaction term between two types of particles 
have two important contributions to the relative ratio of 
energy densities of these particles. The first one has a 
classical interpretation and it is due to the cosmological 
evolution of the different energy densities. This evolu- 
tion will be determined the dynamics which includes the 
potential V and the interaction term Li nt and is calcu- 
lated by solving the equation of motion of the lagrangian 
in a FRW metric. The second consequence of the inter- 
action term L int is to allow for a particle decay, which 
is due to quantum physics. Both effects are relevant for 
determining the ratio of energy densities. 



A. Interacting Dark Energy with Fermions 

The interaction between dark energy and fermions i\) 
can be achieved by taking a Yukawa type interaction [8] 



Lint = -f{<j>)lftP 



(5) 



with a function / that depends on the scalar field <fi. This 
term gives an effective field dependent mass to ip 



M = /(</>)■ 



(6) 



In principle the function / is an arbitrary semipositive 
definite function of (f>, f(<f>) > 0. However, it is useful 
to assume that it is a monotonic function, otherwise not 
much can be said about /. 

The evolution of the energy density of the fermions 
Pij, can be determined using the Fermi-Dirac distribution 
and taking into account that the mass M is a function of 
the scalar field <fi, the evolution of ip and 4> are given by 
[8] 

p$ + 3H p^(l + w^) = (l-3w^,)p^ — (7) 
4> + 3H4> + V = -(l-3«)^j (8) 



with V(4>) the potential for the scalar field 4> in the ab- 
sence of interaction and /' = df/d(j). The solution to 
eq.(7) is 



Pip = p^ a 



-3(l+uty) 



1— 3w^, / j 



(9) 



where from now on the subscript o is present time. In this 
case the fermion fluid no longer redshifts as a^^ 1+w ^ 
since the evolution of f{<j>) also contributes to the red- 
shift. Eq.(9) reduces to p^ = p^ a~ 3 f/f when w^p = 
which was used in [8]. The quantity fty is independent 
on <j>. We see that for radiation w$ = 1/3 the interaction 
term is 5 = (1 — 3w^)p^f / f = and p,p = p^p ~ a~ 4 , 
consistent with having a fermion field with a constant 
(vanishing) mass. From now on we will set = for 
simplicity but the generalization is straight forward. 

Eqs.(7) and (8) can we rewritten in terms of p^, defined 
in eq.(2), and p^ as 



Ptj, + 3Hp^(l + 



f 

Pi>j = S(t), 



(10) 



p <t> + 3Hp <t> {l + w (j ,) = - Pi , J - = -8{t) (11) 

where we have used / = ftp. The interaction term is 
defined by 5 = p^f/f is a function of time and we could 
solve eqs.(10) without making reference to the scalar field 
</>. We can write eqs.(10) and (11) as 



p0 = SHp^l + Wfeff) 
p^j = -3Hp^(l + w^ e ff) 

with the effective equation states given by 



w i>eff = w i> - 

W 4>eff = W <P + 



3Hp^ 
8 

3Hp <j> 



f'4> 

3Hf 

Pip f'4> 

P<p ZHf 



(12) 

(13) 
(14) 



For 5 > we have w e f / > w^ and the fluid p^ will dilute 
faster then without the interaction term while p^ will 
dilute slower since w^, e ff < w$. 

The complete evolution of p^, and p^ depends on the 
effective equation of state parameters defined in eqs.(13). 
Which fluid dominates at late time will depend on which 
effective equation of state is smaller. The difference in 
eqs.(13) is 



Aw, 



eff = W^ e ff - Wfaff 



Aw 



with 



+ Pi> 



3Hf 

While the sum gives the cosntraint 



(15) 



(16) 



QipW^ e ff + fl^Wfeff = fi^Uty + f^W0. (17) 
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Clearly the relevant quantity to determine the relative 
growth is given by T and if T > Aw we have Aw e f / < 
and p^ will dominate the universe at late times while for 
T < Aw we have Aw e /f > and p^ will prevail. 



B. Effective Potential Vr 

Form eq.(8) we see that the effective (total) potential 
Vt and its derivative w.r.t <j> are given by 



V T = V + p,p = V + p^p 



fo 



V T = V + P ^ f - = V + p^t 

J Jo 



(18) 
(19) 



where we have used eqs.(9) in the last equality of 
eqs.(18)and (19). The mass of <f) is given by 

m 2 = V^ = V" + Pi , f — = V' + pj, f —. (20) 

/ Jo 

Using eq.(19) we can rewrite eq.(8) as 

4> + 3H4> + Vt = 0. (21) 

While the derivative of the effective potential V T is 
negative the field <j> will evolve to larger values since 
4> > 0. However, a minimum of Vr can be reached if 
V T — which requires /' > 0, since V < by hypoth- 
esis, with V = —p^f'/f- Taking the time derivative 
V T = V T '<j) + d t V T = one obtains [10] 



3HV _ 3Hp^f _ 3Hpjf 



m 2 f 



m 2 f 



(22) 



Notice that even in the case with V T — we have <j> > 0, 
i.e. the <f> grows with time even at the minimum of Vr, 
and we also have / = f'<f> > 0. The mass m given by 
eq.(20) becomes [10] 



m 2 = V" + pi> J — = Pii, J — I 1 



f-*7( ,+ ^l :> 



with T = V"V/V' 2 (T m > 1 if the field 4> is tracking [3]) 
and we used V = —p^f'/f- Approximating V ~ we 
get 



rrf f_( f 2 - 0+ 

h 2 ~ 3 ^ / [ 1 + ff» Tm nt 



(24) 



Solving eq.(21) with V T = (i.e. the usual slow roll con- 
dition <j) <C 3Hcj) is no longer satisfied) gives <f> + 3Hcf> = 
with a solution <ft ~ a ~ 3 - Since p^, redshifts as a~ 3 we 
see from eq.(22) that 



3Hf _ f y/ptil+w*) 



Pi, m 2 f f 



Pi: 



(25) 



with k a constant and we have used from eq.(3) <j) 2 
p^(l + w^) and p^p = p^f /f. By taking 
p<l,/3H 2 ,Q.^ = p^/3H 2 , eq.(25) with M = f gives 



1 M 2 



fc 2 H 2 fi (1 + W4> ) ' 
Dividing eq.(24) by eq.(26) we get 



(26) 



m 
M 2 



+ n*f" f, , f 2 r Q A (97) 
k 2 / \ 1 + ff" Lm nJ- (27) 



At present time with 0^, o = 0.7, Q$ = 0.3, w<j, = —0.9 
we have f#A^(l + ^ 1.3 and M 2 /H 2 k 2 ~ 5. 



IV. DECAY WIDTH T 

The conditions for a particle to decay is that its lifetime 
t = 1/r, where T is the decay width, is smaller than the 
life of the universe given by t oc l/H, 



(28) 



and that the energy condition (in the rest frame of the 
decaying particle) 



m 



2 > £,M 2 



(29) 



must be satisfied, where m is the mass of the decaying 
particle and the sum is over the product particles with 
mass Mi. 

We assume an interaction term L tnt — — f(<t>)i>i> as 
in eq.(5) with / not necessarily a positive power law 
function of <j>, as for example h cx e b< ^ or h oc 1/0 2 . If 
we expand f in a Taylor series around a time t = t e , 
<(> e = <f>(t e ), we have /'= f e + f'\ e 6<f> + (l/2)/"| e <S</> 2 + 
with /' = df/d(j) and f e — f(<fi e ), then the interaction 
term gives an effective coupling 



-Lint = f{(j>)i>i> 

= felftP + fe &<t> W> + 



(30) 



r:54> 2 w+ 



between two fermions and q quantum scalar fields <fi, with 
q = 1, 2... . The first term in eq.(30) gives the mass of the 
fermion while the second gives an interaction between one 
scalar field and two fermions. This term is responsible 
for the scalar decay. At the time t e the fermion field gets 
a mass given by M — f e . Since the potential of <f> might 
be a run away potential, i.e. the minimum is at <f> — * oo, 
the expansion point <ft e is a function of time and the mass 
M and couplings / e , f' e are also functions of time. 

If / = h(f>, with h constant, the decay rate of <j> into 
two fermions is given by T = h 2 m/8n and generalizing it 
to an arbitrary function f{<f>) we obtain 



r = 



f ,2 m 
8tt 



(31) 
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The dark energy scalar field will decay if Y/H > 1 and 
m 2 > 2M 2 which implies that 



/' 4 m 2 



3/' 4 / V" 



H 2 (8tt) 2 H 2 (8tt) 2 \3H 2 



or 



2M 2 



1 / y" 

2 1/ 



- ( "77 ' P+Ja ) ' [ 



Sl+^r ) > 1 (32) 



(33) 



with m 2 = V" + p^f'/f as defined in eq.(20). 



V. GENERIC CONDITIONS 

The interaction between dark energy and p^ must be 
such that it allows the universe to accelerate recently. 
Therefore the dark energy particle must not decay before 
present day. However, the dark energy may decay in the 
future and the universe would then stop accelerating. We 
have therefore different sets of conditions: for t smaller 
than t , present day t = t Q and for times t > t Q . 



A. Conditions for acceleration 

For the universe to accelerate in the absence of any 
interaction the slow roll conditions in eq.(4) must be sat- 
isfied. However, once the interaction is turned on the 
slow roll conditions change. 

For generality purpose we take the energy density of 
the universe as p-^ + p^, plus p n im, with p n i m a non in- 
teracting matter with w n im = 0. The Hubble parameter 
is then 



3iP 



1 



P^> + Pnim = -z4> + V{4>) + P^> — f 1" P 



fo 



1 > 



+ V T +Pn 



(34) 



Using eqs.(10) and (11) it is easy to see that the universe 
accelerates if 



H 2 + H = ~ (p m + p,p+ p<p + 3p^) 
-i( Pm + 20 2 + V T -3v) 
-i(p m + 20 2 + p^-2u) 



is positive. A positive acceleration gives the constraint 



^ Pnim + P4, + < 
2 V 



< 1. 



(36) 



Clearly V must be larger than (p m + p t p)/2 and p^ > 
(Pm + Pip)/2- Eq.(36) gives the usual constraint il^w,/, < 
— 1/3, even though p^ is also a function of <j>. So as long 



acceleration. For recent times > p^p but well before 
that we will have the dark energy as subdominant energy 
density p^ < p m = p nim + p^ . 

We would like to emphasize that the " naive" conditions 
on the total scalar potential Vr 



\YlL\ 
U T 



< 1 



(37) 



do not imply an accelerating epoch. The reason is that 
only the potential V and not Vr enters with a negative 
sign (within the brackets) in eq.(35). The proof is simple, 
just take the limiting case V = V = V" = 0, in this case 
V T = P4, and V^/V T = f'/f, V^/V T = f"/f can be both 
much smaller than one but condition in eq.(36) is not 
satisfied and therefore the universe does not accelerate. 



B. Conditions for t < to 

For dark energy not to decay before present day we 
need the constraints in eqs.(32) and (33), Y/H < 1 or 
m 2 < 2M 2 , not to be satisfied simultaneously, therefore 
we need 

11 

H 2 

or 



f 4 m 2 3/ M 
(8ir) 2 H 2 ~ (8tt) 2 



V" f" 
3IP +Q *J 



< 1 (38) 



m 

2M 2 



1 fV" f"\ , 

- \j2+P^ ) < 1 



(39) 



Clearly the interaction term /(</>) and its derivatives play 
the crucial roll. The term V"/3H 2 in cq.(38) is much 
smaller than one by hypothesis so the inequality requires 
3/"/'7/(8tt) 2 < l/JV 



C. Present day t = to 

At present day we know that the universe is dominated 
by dark energy with 



f^0o = 0.73, 



^mo — 0.27 



(40) 



and dark matter fi m includes non interacting dark matter 



(35) 

p n im and interacting dark matter p, 



Wi Pm 



+ 



Pi m . In our case the interacting dark matter is given 
by fermions, p im — p^. For simplicity we could take 

Pm — Pij)! Pnim 0. 

The cosmological data for a constant equation of state 
of a dark energy, assuming no interacting with other flu- 
ids, is [? , DE]' 



-1.04 ±0.06. 



(41) 



as Q,p < 1/3, or Q r , 



The equation of state of sate (given by eq.(3)), is 
Wrj, > —1 but an apparent w app < — 1, as defined in 
fi^, > 2/3, there will be no eq.(56), can be obtained due to the interaction between 



■5 



<p and the fermion fields. Since w app is smaller than -1 we 
take this as an indeication od dark energy interaction. 

At present time p$ > and using the condition for 
acceleration in eq.(36) we need 



^1 <2 



(42) 



where we have taken for simplicity p n i m — and the slow 
roll conditions in eq.(4), which imply that <j> 2 /V -C 1. 



D. Conditions for dark energy decay for t > t 

Dark energy will stop dominating the universe if ei- 
ther p^ redshifts slower than p^, due to the interaction 
terms /(</>), or if the dark energy particle decays into the 
fermion fields. 

Dark energy density will dilute compared to the 
fermion energy density if 



T = 



/ ( P<t> + P± " 
3Hf \ P0 



> Aw. 



(43) 



Dark energy decay will take place iiT/H > 1 and m 2 > 
2M 2 given by eqs.(32) and (33). 



VI. FLUID EVOLUTION 

The complete evolution of p^ and p^ depends on the 
effective equation of state parameters defined in eqs.(13). 
The difference in eqs.(13) is [5] 



W <j>eff 



Aw e ff = W^ ef f 

with Aw = w^ — w,p and 

T = J_ ( P4> + P4> 
~ 3H 



Aw-T 



(44) 



(45) 



If we take the ratio y 
of y w.r.t. time is [5] 



P4>P4> 

P<plPip = f^/f^ the derivative 



y = 3Hy [Aw - T] 



(46) 



The value for y is constraint to < y < oo with y = 
for — and y = oo for p^ = 0. Clearly from eq.(46) 
we see that the evolution of y depends on the sign of 
Aw — T. Clearly the relevant quantity to determine the 
relative growth is given by T and if T > Aw we have 
Aw e ff < and will dominate the universe at late 
times while for T < Aw we have Aw e ff > and p^ 
will prevail. For no interaction (5 = and T = gives 
Aw e ff — Aw > and p^ dominates at late times. If 
T = Aw then w^ e ff — w<f, e ff and the ratio of both 
fluids p^/pcj, will approach a constant value, and if the 
universe is dominated by p^ + p^, i.e. fl^ + Q^, = 1, then 
eq.(17) gives 



W0 < w ef f = w^Vl^ + w<f,(l - fi^) < w^, 



(47) 



i.e. the effective equation of state is constraint between 
W0 and w$. 



1. Non Interaction solution: Aw > T 

If (tity, — w<p) > T then y is positive and y will increase, 
i.e. pip will dilute faster than p^, and we will end up with 
p<p dominating the universe. In this case the interaction 
term <5 is subdominant and the evolution of p^ and p^ is 
the usual one, i.e. oc a^ 3 ( 1+M, «) and p^ oc a^ 3( - 1+w ^\ 



2. Finite solution: Aw = T 

A solution to eq.(46) with y constant and p^ ^ 0, p^ ^ 
is only possible if S/Hp^ is positive and constant since 
T = (5(1 + y)/Hp ( f > = Aw must be constant and positive, 
taking Aw = w^ — w^ > constant. Let us take S = 
CIIp^ with C a positive constant, and from eq (46) for 
y = we get a stable value of y given by [5] 



_ ^0 _ 3(w^ - w ) _ 



(48) 



It is easy to see that the solution y s is stable since from 
eq. (46) the fluctuation Sy = y — y s to first order behaves 
as Sy/Sy = —CHy s (CH is positive by hypothesis) giving 
Sy — > 0. The evolution of p^ is given by [5] 



-3(l+«v)-C 



-3(1+- 



P<p ~ P<t>i\ — ] — P<t>i 

\<H / \ u,j / 

(49) 

where we have used that C = 3(w^ — w^) /(l + y s ) . Since 
C is positive the solution in eq.(49) dilutes faster than 
the non interacting solution p^ oc a^ 1+w ^ and with an 
effective equation of state w^ e ff = w^ + C/3 = (w$ + 
ysW^/il + ys) > w . 

3. Interacting solution: Aw < T 

For p^p to dominate the universe we need y <C 1 at late 
time and the (interaction) term T should dominate over 
Aw. A simple example is when A — 5/p# is constant and 
positive. In this case the evolution of y is given by 



a 

Vi \ — 

ai 



3(w^,-w^) 



-At 







(50) 



W0. 



at late times for any value of w^ 
A as a constant decay rate T. 
I) Constant decay Rate T = A = 8 / p^ 



Now, let us take 



If the mean lifetime of a particle is given by a constant 
r then we can expect it to decay into lighter fields. In 
the case T = 1/t constant the constrain T > 3H will be 
satisfied at some time i ^ r (or equivalently at T/H ~ 1). 

For a constant T the interaction term is 6 = Tp<p and 
the solution to eq.(ll) using (14) is simply 



a 

P4>i I — 



-3(l+uv) 



-r< 



(51) 
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For small t, i.e. for 3H(l + wp) ^> T, one has the 
usual redshift in the absence of any interaction pp oc 

a -3(l+w+) K X j t 2 wmlc for large tj j g 3^ (J + ^) ^ Fj 

one has an exponential decrease pp oc e~ tr . Since T is 
constant, then clearly Tp^/p-p, will not be constant and 
the solution to eq.(ll) will be quite different to that of 
P<t>- 

Let as take the ansatz p,p, = p^ a^ 1+w ^ + g(t)p^ 
taking the time derivative we have 

f>$ = -3H(l + w^)p^ + gp<t> + gp<p (52) 
= -3JT(1 + u^W + (3H Awg - Tg + g) P<f> 

and eq.(52) should be equal to eq.(10), i.e. T = S/pp = 
3HAwg + g — Tg. Therefore, the function g must satisfy 
the equation 



3AwHg + g 
1 + 5 ' 



(53) 



Eq.(53) is valid for any functional form of T and not 
necessarily a constant. However, a simple solution can 
be found when T is constant for small t. The solution is 
g = git giving T ~ (3AwHi + l)gi constant, and using for 
the dominant energy density Hi = 2/3(1 + wp) we have 
T = <7j(l + 2wjf, — w<p)/{l + wp,). The solution to eq.(52) 
is then 



= P^ a- 3{1+W ^ 



Pi> — Pi> 



(54) 



with q = (1 + wp,)/(l + 2w^ — wp) (if we have matter 
decaying into radiation, w<p = and = 1/3, then 
q = 3/5). Eq.(54) shows that the fluid redshifts as usual 
for Tt <C 1, i.e. p^ ~ a -3( 1 + M, V') j where the first term 
in eq.(54) dominates. At t ~ ti(Yti)~ q the second term 
starts to dominate giving p^ oc p^t oc t^ 1 oc a _3 ( 1 + M '«)/ 2 j 
with the fluid decaying into p^ at around Tt ~ 1 when 
p^ starts decreasing exponentially fast. For Tt 3> 1 the 
fluid p^, will dominate the universe redshifting again as 

Pi , - a -3(l+«v). 

We have seen that a constant decay rate T gives an 
exponentially suppressed energy density p^. However, a 
constant decay rate is not realistic. The decay rate T for 
a scalar field decaying into fermions, given in eq.(31), is 
a function of m and /' and both are in general field and 
time dependent so T will not be constant. 



A. Apparent Equation of State 

An interesting result of the interaction between dark 
energy " DE" with other particles is to change the appar- 
ent equation of state of dark energy [6]-[10]. An observer 
that supposes that DE has no interaction sees a different 
evolution of DE as an observer that takes into account for 
the interaction between DE and another fluid. This ef- 
fect allows to have an apparent equation of state w < — 1 
for the "non-interaction" DE [10] even though the true 
equation of state of DE is larger than -1. 



Let as take the energy density p = p^ + p^p = Pde + P4>- 
The energy densities are given by eqs.(ll) and 

(10) and these two fluid interact via the S term. On 
the other hand the energy densities pde and p^, do not 
interact with each other by hypothesis and therefore we 
have p^p = —3H (we have taken = 0) and puE — 
-3Hp DE {\ + w ap ), i.e. 



Pde 



PDEo a 



-3(l+w ap ) 



(55) 



if-- 



'app 



is constant. It was pointed out that the apparent 



equation of state w app can take values smaller than -1 
and it is given by [10] 



W ap = 



P4> \fo 



wp 
1 -a 
_P± 

P4> 



(56) 



P±[ x _[o_ 

./• 



and we have used eq.(9). We see from eq.(56) that for 
P4> < p^p, i.e. f(t) < ,fo(t ) for t < t a , we have x > and 
w ap < wp which allows to have a w ap smaller than -1 in 
the past. However, for / > f with x < and w ap > wp 
for t > t Q , The observational evidence shows that / is at 
present time growing function of <fi, i.e. /' > 0. 



B. Evidence for Dark Energy Decay 

The cosmological observations prefer an equation of 
state w < — 1 for dark energy. In principal a w < — 1 
for a fluid is troublesome since it has instabilities and 
causality problems. However, as seen in section VI A this 
can be an optical effect due to the interaction between 
dark energy with other particles as for example fermions. 
It was shown that for an effective w < — 1, interpreted 
as an interaction between dark energy and another fluid, 
can be a signal for dark energy decay in the future with 
the universe no longer accelerating [5]. 

Lets us expand / around present time t f = f a + f a St 
and we keep only the first order term. Now, from eq.(56) 
x can be written as [5] 



P4> \Jo 



_P±fo 5t = _P±fo 5t 



fo 



P4> f 



(57) 



with pip = p,pf /,f, St = t — t a . For t < t we have x > 
since f > 0. The interaction term 5 = p-pf /f 7 is now [5] 



S = -p^ 



St 



(58) 



where we have taken / = f a , consistent with the approx- 
imation taken for /. Let us take the simple phenomeno- 
logical ansatz for x given as x — {3(a — a) = —[35a with (3 
a constant to be determined by observations. The SNla 
data are in the range 1 > a > 2/5, i.e. for a redshift 
< z < 1.5, and the best fit solution has an average 
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equation of state < w >« —1.1 [1]. Taking the average 



of 



we have 



J 01 Wq 



, da 



f a ° da 

Jai 



w^Log [1 - (5(1 - ai) 
/3(ai - 1) 



(59) 



where we have used a; = —/35a and a a = 1. As an ex- 
ample let as take < w apP >= —1-1, as suggested by the 
observations [1, 2], [9], and wj, — —0.9, tity = 0, ai = 2/5. 
In this case we obtained from eq.(59) the value (5 = 0.56 
[5]. Using the anstaz x — — (58a with 8a = aHSt then 
eq.(58) becomes [5] 

8 = (5aH PM> , L = /3aH and T = ^ ( 1 + — J • (60) 
J 3 \ P4>J 

Using cq.(ll) with the interaction term given in eq.(60) 
we get an energy density [5] 



a -3(l+uv) e -/3(a-l) 



(61) 



which shows that dilutes as a~ s ( 1+w ^ for a <C a a = 1 
and is exponentially suppressed for a > 1. In this 
case 3Aw < T and T — > oo at late times, as seen from 
eq.(60). 



VII. DARK ENERGY DISAPPEARANCE 

The decay rate T = f'm 2 /87T, the mass m and M 
are a field and time dependent quantities. If we take 
/ as a monotonic function, since /' at present time is 
positive (cf. section VI B), then / should be an increasing 
function of cf) for all values of <fi. This implies that the 
mass of the fermion M — f increases at all times. So, 
if the mass of the scaler field m vanishes asymptotically 
then m/M will be smaller than one at some point in the 
future but larger in the past. So, we would expect on 
general grounds that the condition m 2 > 2M 2 would be 
easily satisfied in the past. Therefore we require T/H < 1 
and from eq.(38) we need 



3/"/ 



< 1 



(62) 



where we have taken p^ 3> p<p (fl^, = 1) in the far past. 



A. Dark Energy Dilution 

As we have seen in section VI dark energy dilution, 
P<p ^ Pip at late times, will only take place if T > Aw 
with 



T 



/ ( P<f> + P<P 



3Hf 



P<t> 



(63) 



Let us take f/fH — (a/ f)(df /da) with the ansatz / = 
f Q a q . We then have f/fH = q, T = g/3 and the in- 
teraction term 8 — qp^H . For q > q a = 3Aw we have 



w 4>eff < w 4>eff a n d p^> will dominate while for q < q Q 
then will end up dominating the universe. Further- 
more, if / grows faster than power law then T — > oo 
and p^ will dominate at late times, this is the case dis- 
cussed in section VI B as sen from eq.(60). For / growing 
slower than power law we have T — ► and p^ will end 
up dominating the universe. 

The limiting case is for p^ / p^ constant and a solution 
is obtained for T = g/3 = Aw. This case is equivalent as 
in section VI 2 with p^ <-> p^ and q = C. From eq.(48) 
we have 



P± 

P<t> 



fl^ 
fid, 



3Az 



(64) 



(see discusion belown eq.(48)). 

Now, we consider again the interaction term f(<j>) as a 
function of a scalar field <f> with / = f'<j>. The evolution 
of 4> is determined by eq.(7) and the sign of V^, defined 
in eq.(18), sets if <f) grows or decreases with time. We will 
discuss the cases Vj. > 0, < and = separately. 



1. Case: Positive Derivative V? > 

First, let us consider the possibility that > which 
implies from eq.(21) that <j) < 0. In this case since V < 
then /' > and / = ftp is negative giving a negative T 
(cf. eq.(63)) and p^ will dilute slower than (even 
slower than without the interaction) . 



2. Case: Negative Derivative V? < 

For Vj. < we have \V'\ > p^f / f and from eq.(8) we 
can approximate 

4> ~ -V/3H. (65) 
Taking, for simplicity pd, + p^p = 3H 2 in eq.(63) we get 



P4>f 



f V 

-±-—>Aw. 
f 



(66) 



The condition V? < implies /'// < \V'\/p^ and to- 
gether with the slow roll conditions p^ > V > \V'\ and 
eq.(66) we get the constraint 



n;>7 >3A *" 



which sets an upper value 
1 



fl-4, < 



< 



1 + /'// 1 + 3Aw ' 



(67) 



(68) 



If we take for example w$ — and Wd, ~ —1 we have 
Aw = 1 and fl^ < 1/4. Therefore, we conclude that in 
the case V T < it is not possible for p^ to dominate the 
universe. 
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3. Case: V? = 

As we have seen, while the derivative of the effective 
potential is negative the field 4> evolves to larger values 
since <ft > 0. However, a minimum of Vr can be reached if 
= 0, which requires /' > since V < by hypothesis, 



and V = —pi/jf'/f- Taking the time derivative V? 
one obtains with 

• _ 3HV _ -AHp^f' 







m 2 f 



(69) 



(c.f. eq.(22)) which is still positive, i.e. even tough V? — 
the field 4> still grows with time. The reason being that 
ffy is a function of time (through a(t)). Calculating T = 
Htf'/pjf > Aw and using eq.(69) and 1 > \V'\/V = 
Pipf'/ fV we get the constraint 



Aw 



which sets an upper value 



£0 



< 



l + (m 2 /'3H 2 )Aw' 



(70) 



(71) 



If we take fi^ + fi^ = 1 then eq.(72) gives un upper limit 



0,0, 



1 we 



fl^p < 1/(1+ \ f" / f\) and in the case 
require \ f"/f \ -» 0. 

The constrain in eq.(32) gives V" > H 2 /f 4 (8n) 2 and 
using iJ 2 > V/3 and V > V" we obtain 



i 

24^ 



(73) 



(74) 



while for condition in eq. (33) we have 

Since V" — > by hypothesis, eq.(74) implies that / — > 
at late times, i.e. / is a decreasing function at large £, 
since M — f > 0. Eq.(74) can only be satisfied if /' 
is negative, i.e. / is a decreasing function. However, 
at present time the dark energy interacting must have 
fg > if w < — 1, as seen in section (VI B). So, unless / 
increases at present time and than it starts decreasing at 
a later stage (i.e. it is no longer a monotonic function of 
(f>) the conditions for a dark energy decay into fermions 
with V" > Pipf" ' I f cannot take place. However, / will 
no longer be a monotonic function of <\>. 



2. Case: m 2 ~ py,/"//, i.e. p^f'/f > V" . 



Comparing <j> form eqs.(65) (labeled 4>i) and (69) (labeled 
</>2) we see that the ratio 4>i/4>2 — m 2 /9H 2 is exactly that 
in the denominator of eqs.(68) and (71). For 4>\/4>2 — 
m 2 /9H 2 <C 1 the evolution of <f> is faster in the case 
V T = than for V T < 0. Now, if m 2 /3H 2 > 1 and taking 
Aw = 1 we have O^, < 1/2. However, if m 2 /3H 2 -C 

« 1 



4> 



1 than then the constraint in cq.(71) becomes 
allowing for the fermion fields to dominate at late times. 

It is easy to see from eq.(70) that in the limit Sl^ — > 
we need /'// -> and m 2 /3H 2 -> 0. From eq.(23) the 
condition m 2 /3H 2 — > implies that /"// — » since 
= 1 - ^0 = 1 and V"/V < 1 by hypothesis. 



B. Dark Energy Decay 

The conditions for dark energy to decay into fermions 
given in eqs.(32) and (33) are constraints on the dark 
energy mass m 2 = V" + p^f"/f- Let us consider the 
two possible cases V" > p^\f"/f \ and V" < p^\f" / f\ 
separately. 



1. Case: m 2 ~ V" , i.e. V" > Pj,\f"/f\. 

Using the slow roll hypothesis on V given in eqs.(4) we 
can set the upper limit p<$, > V > V" > p^\f"/f \ which 
gives the constraint 



fi,i, t 



(72) 



Let us now take the case m 2 ~ p^f"/ f > V", which 
implies that 



n ^ V " PtP Pip 



(75) 



where we used eq.(9). Eq.(75) shows that at late times 
(with a — ► oo ) we have V" / f" — ► 0, i.e. V" must decrease 
faster than /". Using m 2 ~ P4>f"/f and condition in 
eq.(33) we get in this case the constraint 



0<J- < P± 
~ I" 2/ 



(76) 



where we have used again eq.(9). Once again we see that 
at a — > oo that ///" — > at late times. Finally, the 
constraint in eq.(32), with m 2 ~ p^f" / f and 3H 2 > p^, 
gives 

rr 



/ 



> 



3(8^) 2 



(77) 



Since from eq.(76) ///" must vanish at late times, 
eq.(77) is satisfied as long as /' 4 does not go to zero 
faster than ///". 



C. Examples 

To have dark energy decaying the conditions V" / f" — > 
0, ///" and f"f' 4 /f > 1/3(8tt) 2 must be satisfied. 
A simple example is 



/ = f e a « 



(78) 
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with a a positive constant. Since <f> grows with time, 
so docs /,/',/" and the constraint V" / f" <C 1 and 
f"f 4 /f > f 4 > 1/3(8tt) 2 are satisfied trivially and the 
condition 

■L = 2a + Aa 2 4> 2 > 1 (79) 

is also satisfied for cf) — > oo. Clearly, this potential will 
not work for dark energy dilution. 

Dark energy diluting needs /'// — > and /"// — > 
and an example is 

/ = M a (80) 

with a > 2. In this case /'// - 1/0 and /"// - 1/0 2 
and grows with time. 

D. Summary 

To conclude, we have seen that if the conditions on 
the interaction term /'// — ► and /"// — ► are sat- 
isfied then we can have a dark energy diluting faster 
than the energy density of fermions. In this case we 
will have m 2 /H 2 <C 1. On the other hand, for dark 
energy to decay into fermion fields we require that dom- 
inant part of the scalar mass is given by the interac- 
tion term m 2 ~ Pipf"/f V" . In this case, the condi- 
tions to be satisfied are given by eqs.(75), (76) and (77), 
i.e.V"/f" - 0, ///" - and J"/' 4 // > 1/3(8tt) 2 . 

The condition for dark energy redshifting faster than 
fermions /"// <C 1 and the condition for dark energy 
decay into fermions ///" <C 1 cannot be simultaneously 
meet. Clearly form eq.(36) if the universe is dominated 
by than it will not accelerate. We have also seen 
simple examples of the interaction term / which require 
no fine tuning of any parameter. 

VIII. CONCLUSIONS 

A dark energy component is responsible for the present 
stage of acceleration of our universe. If no fine tuning is 
assumed on the dark energy potential it is easy to see that 
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since it dilutes slower than the other fluids, e.g. matter, 
then it will end up dominating the universe at late times 
and the universe will not stop this stage of acceleration 
anymore. 

In this letter we have studied the possibility that the 
universe will stop accelerating and that dark energy de- 
cays into fermion fields. The interaction between dark 
energy and fermions is given through a fields dependent 
mass M = /(</>) of the fermion fields. 

Furthermore, the fact that the equation of state of dark 
energy seems to be smaller than minus one as suggested 
by the cosmological data can be an indication that dark 
energy does indeed interact with other fluids. We take 
this fluid to be fermions. Using the observational data we 
can determine present day value of the interacting func- 
tion f(<f>) and its derivative and we can then extrapolate 
the result into the future. 

The interaction term /((/>) plays an important roll in 
the evolution of the dark energy. We determine the con- 
ditions under which a dark energy can dilute faster than 
the fermion fields or it can decay into these fields. We ob- 
tained that only in the case where the mass of the dark 
energy, given by m 2 = V" + pi/jf"/f, is dominated by 
the interaction term can the dark energy decay. While 
for dark energy diluting faster then the fermion fields the 
conditions /'// 

rightarrowO and /"// — > must be satisfied. 

We have shown that naive slow roll conditions on the 
effective potential Vr do not imply an accelerating epoch. 
The condition needed is that the scalar potential V dom- 
inates the universe (cf. eq.(36)) 

We have seen, therefore, that it is indeed possible to 
live now in an accelerating epoch dominated by the dark 
energy and without introducing any fine tuning parame- 
ters the dark energy can cither dilute faster or decaying 
into fermions in the future. The acceleration of the uni- 
verse will then cease. 
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